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ABSTRACT 


In  this  study,  the  author  uses  a  MATLAB  simulation  to 
develop  and  test  a  generalization  of  the  traditional  Random 
Search  model  which  allows  both  the  searcher  and  target  to 
move  and  to  be  in  different,  but  overlapping,  areas.  Also 
the  best  evasion  speed  for  a  randomly  moving  target  against 
a  Systematic  Search  is  studied. 
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EXECUTIVE  SUMMARY 


In  this  study,  the  author  uses  a  MATLAB  simulation  to 
develop  and  test  a  generalization  of  the  traditional  Random 
Search  model  which  allows  both  the  searcher  and  target  to 
move  and  to  be  in  different,  but  overlapping,  areas.  Also 
the  best  evasion  speed  for  a  randomly  moving  target  against 
a  Systematic  Search  is  studied. 

The  new  generalized  Random  Search  formula  is, 

Ft  (i t )  =  aFTslow(t )  +  (1  -  a)Frfml  (7)  ,  where 
a  =  exp(-^Vo  , 

E  A,U 

4AARV  ' 

Ft,, At)  =  ( A,  /  4)(  1  -  exp(-2 RVll A,))  , 

FTfas, (0  =  1- exp(-2i? V (Ast  IAt)tl  As ) ,  and 
V  =  mean  relative  speed  between  searcher  and  target. 

With  regard  to  the  best  evasion  speed  against  a 
Systematic  searcher,  extensive  simulation  suggests  that  a 
5%-20%  of  searcher  speed  is  the  best  speed  for  the  target. 
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I. 


INTRODUCTION 


A.  BACKGROUND 

Random  Search  is  a  popular  model  for  area  search 
because  it  is  mathematically  simple  and  it  provides  a 
conservative,  lower  bound  on  the  probability  of  detecting  a 
stationary  target  with  Systematic  Search.  However,  the 
Random  Search  model  also  has  significant  limitations.  In 
particular,  it  assumes  that  the  searcher  and  target  are 
contained  in  the  same  area  and  that  the  target  is  stationary 
We  address  these  model  limitations  in  this  thesis. 

B.  RESEARCH  QUESTIONS 

This  thesis  will  address  the  following  two  questions: 

1.  How  can  the  Random  Search  model  be  generalized  for 
situations  where  the  searcher  and  target  areas  are  not 
coincident  and  both  the  searcher  and  target  are  moving? 

2.  What  is  the  best  speed  for  a  randomly  moving  target 
to  evade  a  searcher  conducting  a  Systematic  Search? 

C.  THESIS  ORGANIZATION 

The  thesis  is  organized  in  the  following  manner: 

•  Chapter  II  reviews  the  search  models  of  Exhaustive 
and  Random  Search. 

•  Chapter  III  presents  a  MATLAB  simulation  of  Random 
Search . 
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•  Chapter  IV  develops  the  extended  Random  Search 
model  where  searcher  and  target  areas  are  not 
coincident  and  both  the  searcher  and  target  are 
moving . 

•  Chapter  V  investigates  the  optimal  speed  for  a 
target  evading  a  Systematic  Search. 

•  Chapter  VI  summarizes  the  results  of  simulations 
and  recommends  future  studies. 
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II.  OVERVIEW  OF  RANDOM  SEARCH  THEORY1 


A.  EXHAUSTIVE  SEARCH 


1.  Characteristics  of  Exhaustive  Search 


Exhaustive  Search  is  an  idealized  search  model  which 
assumes  a  stationary  target  and  a  'perfect  search''  ,  meaning 
no  search  overlap,  no  search  effort  placed  outside  the 
target  area,  and  all  of  target  area  is  completely  covered  by 
the  search  sensor.  The  searcher  has  a  'cookie-cutter  sensor' 
with  range  R  ,  sometimes  called  a  definite  range  law  sensor. 
Such  a  sensor  always  detects  a  target  within  a  specified 
range  R  ,  and  never  detects  targets  beyond  that  range.  That 
is  : 


JP(detection)=P/  = 


[  1,  range  <  R 
[0,  range  >  R  . 


(1) 


Because  of  its  optimistic  area  coverage  assumptions. 
Exhaustive  Search  is  generally  assumed  to  provide  an  upper 
bound  on  the  performance  of  realistic  search. 


2.  Exhaustive  Search  for  a  Uniformly  Distributed 
Target 

Assume  that  a  searcher  conducts  Exhaustive  Search  over 
the  search  area  A  ( inn2 )  with  a  cookie-cutter  sensor  of  range 
R  (  nm  )  ,  speed  V  (  nm/hour  )  ,  and  sweep  rate  2  VR  (  mn/hour  )  .  The 


1  The  references  for  this  section  are  the  unpublished  lecture  notes 
of  Professor  James  N.  Eagle  and  Professor  Alan  R.  Washburn's  'Search  and 
Detection'  book  (2002). 
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time  of  detection  in  the  random  variable  T  (  hour  )  ,  which  is 
uniformly  distributed  between  time  0  and  AK2VK)  .  That  is, 
T~  U  [  0,  A/(2VR)  ]  . 

The  probability  density  function  and  cumulative 
distribution  function  for  random  variable  T  are: 


f(2VR)/A,  te[0,A/(2VR)] 
[O,  otherwise 


fO,  t  <  0 


FT(t)  =  \ 


(2 VRt)IA,  te[0,AI(2VR)\ 


[\,t>AI(2VR) 


(2) 


The  mean  time  to  detection  is  given  by  the  expected 
value  of  T  : 

E(T)  =  A/(4VR)  .  (3) 


Figure  1.  Density  and  CDF  for  the  Time  of  Initial 

Detection  T  with  Exhaustive  Search. 


B.  CONTINUOUS  SEARCH 

1.  Model  Assumption  and  Definitions 

•  y{t)  =  detection  rate  at  time  t  (units:  1/time)  . 

Detection  rate  has  two  equivalent  interpretations: 

4 


(4) 


i.  /’(detection  in  [t,t  +  At])  =  y(t)At  +  o(At) 

~  y{t)At,  for  small  At 

Note:  o(At)  is  a  function  of  At  which  goes  to  0 
faster  than  linearly.  That  is: 


lim 

A  t — >0 


o(At) 

At 


=  0 


rl+At 

ii.  E(#  detections  in  [t,  t  +  At])  =  y(s)ds 

Jl  •  (5) 

«  y(t)At,  for  small  At 

T  is  the  random  time  of  initial  detection,  and  we 
define  FT{t)  =  P(T  <  t),  and  GT(t)  =  P(T  >  t)  =  l-FT(t)  .  (6) 


Detection  events  in  non-overlapping  time  intervals 
are  probabilistically  independent.  So, 

P{N(t  +  At)  =  0}  =  P{N(t)  =  0,  N(t  +  At)  -  N(t)  =  0} 

=  P{iV(t)  =  0}-f5{iV(t  +  At)-iV(t)  =  0}  (7) 

where,  N(t)  be  the  random  number  of  event  occur  during  (0,  /]  . 


2.  Derivation  of  GT(t) ,  FT(t)  and  E(T) 


Since  the  events  of  detection  (and  non-detection)  are 
independent  in  non-overlapping  time  intervals, 

GT(t  +  At)  =  GT{t)(\-[y{t)At  +  o(At)])  .  (8) 

A  B 


Where : 


A  =  P  (non-detection  from  time 
B=P  (non-detection  from  time 


to  0  to  time  t) 
t  to  time  t  +  At)  . 


Re-arranging  terms. 
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Fr  (t)  =  1  -  Gt  (0  =  1-  exp(-  [  y (s)ds) 

J  0 

fr (t)  =  ^-FT(t)  =  y(t) exp(-f  y(s)ds)  .  (10) 

at  J0 

poo  poo  fit 

=  J  ntfT{t)dt=  I  t  K0exp(-l  j{s)ds)dt 

Jl= 0  Jl= o  Jo 


If  J  fT(t)dt  =  \  ,  then  T  is  a  proper  random  variable, 

eventually  detection  is  certain,  and  a  slightly  simpler 
expression  is  possible.  For  any  proper  non-negative  random 
variable  T  : 


fi  CO  poo  fit 

E{T)  =  j  (l-FT(t))dt  =  j  exp(-j  y(s)ds)dt  .  (11) 

The  proof  is  sketched  below  in  Figure  2. 
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Figure  2.  E(X)  for  Proper,  Non-Negative  Random 

Variables  X 

C .  RANDOM  SEARCH 

1.  Characteristics  of  Random  Search 

Assume  a  search  where: 

a.  Searcher  has  a  cookie-cutter  sensor. 

b.  Each  small  segment  of  the  searcher's  track  is 
randomly  and  uniformly  distributed  over  the  search 
area . 

c.  No  search  effort  falls  outside  the  search  area. 

d.  Target  position  is  fixed  in  search  area. 

During  any  time  interval  At  ,  the  searched  area  is 
2 RV At  .  Since  this  area  is  uniformly  distributed  over  the 
search  area  A  ,  the  probability  of  its  covering  the  target 
is  : 

IRVAt  n.,  x  .  ,  .  r 

- =  P( detection  during  [/,  t  +  A/])  . 

A 
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Thus,  we  have  a  constant  detection  rate  with  y  =  {2RV)!A. 


Therefore : 

Ft  (t)  =  1  -  exp(-2i?  Vt  /  A) 

fT(t)  =  (2RV / A)exp(-2RVt/ A)  .  (12) 

E(T)  =  A/(2RV) 

2 .  Dynamic  Enhancement  of  Random  Search 

Now  assume  the  target  with  speed  U  and  the  searcher 
with  speed  V  move  randomly  and  independently  over  area  A  . 
The  searcher  sweep  width  is  2  R  ,  and  T  is  the  random  time 
of  initial  detection.  As  before,  we  wish  to  compute 
Fr(t)  =  P(T<t )  ,  and  for  U « V  ,  we  would  expect  that 

Fr(t)  ~ \-e~2RU/A  .  And  more  generally,  we  will  look  for  a  speed 

F>max(F,{7)  such  that  Fr(t)& \-e~2RUIA  . 


Figure  3.  Dynamic  Enhancement  of  Random  Search 
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Referring  to  Figure  3,  V(9)  is  the  relative  speed 

between  searcher  and  target  when  9  is  the  difference  in 
their  courses.  By  the  Law  of  Cosines,  V  (9)  =  y] V 2  +  U2  -  2  UV  cos  (9)  . 
We  now  assume  that  9  is  uniformly  distributed  between  0  and 
2 n  ,  so  the  average  V (9)  is: 

V  =  fs:V(9)f(9)d9 

=  r  JV2  +U 2  -2UVcos(9)( — )d9 

2  n  (13) 

=  -  JJ  q  ^V2  +U2  -2UVcos(9)d 9 
>  max(F,  U )  . 

This  equation  suggests  that  in  Random  Search  with  dynamic 
enhancement,  the  searcher  speed  is  effectively  "enhanced"  by 
the  target  speed.  The  author  will  confirm  this  by  MATLAB 
simulation  in  Chapter  III. 
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As  shown  above  in  Figure  4,  we  now  consider  the 
situation  where  the  search  area  is  not  completely  overlapped 
with  the  target  area.  Two  types  of  target  behavior  were 
assumed  in  order  to  get  a  simplified  mathematical  expression. 

a.  Slow  Target 

Assumptions : 

-  Target  initial  position  is  uniformly  distributed 
over  At  . 

-  A  target  starting  inside  (outside)  A s  will  remain 
so  for  the  entire  search  time. 

In  this  case  FT(t )  is: 

Fr(t)  =  /’(detection  by  time  1 1  target  start  in  Ast)P( target  starts  in  Asl)  + 

/’(detection  by  time  1 1  target  starts  outside  Ast)P( target  starts  outside  in  Ast) 

=  /’(detection  by  time  1 1  target  starts  in  Ast  (/(target  starts  in  Ast)  (14) 

=(—)(! -exp(-2  RVt/As)). 

A 

b.  Fast  Target 

Assumption : 

-  Target  spends  ( Ast/At )  of  the  search  time  inside 

A  • 

By  time  t  ,  this  target  has  been  available  for 
detection  for  (Ast/At)t  time  units.  Therefore, 

Fr(t)  =  l-exp(-2PV(Ast  /  At)t  /  As)  .  (15) 
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In  light  of  these  assumptions,  eventual  detection 
for  the  fast  target  is  assured.  In  contrast,  the  upper  bound 


of  FT(t )  for  the 
will  attempt  to 
to  allow  FT(t )  to 


slow  target  is  (Ast/At)  .  In  this  thesis,  we 
generalize  the  fast  and  slow  target  models 
be  estimated  for  intermediate  target  speeds 
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III.  SIMULATION  OF  RANDOM  SEARCH 


A.  DESCRIPTION  OF  RANDOM  SEARCH  MODEL 

1.  Characteristics  of  the  Searcher 

The  searcher  is  assumed  to  have  a  "cookie-cutter" 
sensor  of  radius  R  .  In  addition,  the  searcher  is  not 
allowed  to  search  for  a  target  outside  the  search  area. 
Therefore  the  searcher' s  allowable  position  is  limited  by 
sensor  radius  R  .  For  example,  if  the  search  region's  X- 
axis  length  is  50nm  and  the  Y  -axis  length  is  also  50nm, 
then  the  overall  search  area  (  As  )  is  50x50nm2  .  But  the 
searcher's  actual  moving  area  (A’)  is  (50-i?)x(50-i?)nm2  . 

The  searcher's  initial  position  is  uniformly 
distributed  inside  of  the  search  region  A'  .  After  that,  the 
searcher  chooses  his  course  randomly,  independent  of  the 
target's  movement.  The  course  change  event  is  determined  by 
Poisson  process  with  rate  A  .  It  is  assumed  that  the  speed 
of  the  searcher  is  always  faster  than  that  of  the  target. 
This  is  allowed  because  the  searcher  and  target  roles  can  be 
reversed  in  the  simulation,  resulting  in  the  same 
probability  of  detection. 

2 .  Characteristics  of  the  Target 

The  target's  initial  position  is  uniformly  distributed 
over  the  search  region  A  .  The  logic  of  the  target  movement 
is  the  same  as  that  of  the  searcher,  and  the  target  has  its 
own,  independent  Poisson  process  with  course  change  rate  At  . 
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B. 


COMPUTER  ALGORITHM 


1.  Inputs  [Units] 

•  Number  of  simulation  replications,  N  = 500. 

•  Maximum  simulation  time,  ?max=150  [hour]. 

•  The  length  of  search  area  of  X  direction,  /v  =150 
[nm]  . 

•  The  length  of  search  area  of  Y  direction,  /v  =150 
[nm]  . 

•  Searcher  speed,  V =200  [nm/hour] . 

•  Target  speed,  U  [nm/hour] . 

•  Searcher's  detection  range,  R=  2  [nm]  . 

•  Searcher's  course  change  rate,  As  [1/hour]. 

•  Target's  course  change  rate,  At  [1/hour]. 

•  The  unit  time  of  simulation.  At  [hour] . 

•  The  size  of  Search  area,  As=  1502  [  nm2  ]  . 

2 .  Functioning  of  the  Program 

When  a  new  replication  begins,  the  initial  positions  of 
the  searcher  and  the  target  are  chosen  from  a  Uniform 

Distribution  over  the  search  area.  The  only  difference 

between  position  selection  logic  of  the  searcher  and  the 
target  is  caused  by  the  searcher' s  sensor  range  (  R  )  .  In 
order  to  prevent  over-searching,  the  searcher' s  initial 

position  should  be  limited  to  search  area  (A')  . 

The  initial  course  is  also  drawn  from  a  Uniform 

Distribution  between  0  and  2 n  .  The  subsequent  course 

changes  for  searcher  and  target  occur  according  to  Poisson 
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Processes  with  rate  As  and  At  respectively.  In  particular, 
at  each  time  step  At  ,  if  Uniform  _Random(0,Y)  <  As  At  ,  then  a  new 
random  course  Cs  =Uniform  _Random(0,Y)x27r  is  selected  for  the 
searcher.  The  course  changes  for  the  target  are  determined 
in  the  same  manner  by  using  At  . 

When  the  searcher  or  target  encounters  an  area  boundary, 
then  a  random  reflection  occurs.  After  the  reflection,  the 
new  course  in  radians  is  Uniform _Random(C±- .5,  C±+  .5)  ,  where  C± 
is  the  perpendicular  course  from  the  reflection  boundary. 
The  parameter  .5  was  determined  experimentally  to  prevent 
both  "corner  capture"  and  too  many  near  "perpendicular 
reflections."  (see  Figure  6) 
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Figure  6.  Random  Movement  Behavior  at  the  Edge 
Depending  on  Different  Parameters 

At  each  time  step,  we  store  the  position  of  the 
searcher  and  the  target,  and  then  measure  the  distance 
between  them. 

In  order  to  closely  approximate  a  continuous  simulation. 
At  should  be  small.  However,  too  small  a  At  means  too  many 
calculations,  which  in  turn  require  an  excessive  time  to 
simulate.  How  small  a  value  of  At  is  sufficient  to  produce 
an  accurate  result? 


16 


V=200,  U=100 


Figure  7 .  Fr(t)  for  Various  Values  of  At 

R 

Figure  7  above  shows  that  A t ~  is  an  appropriate 

value  to  use.  Thus  the  searcher  moves  half  the  distance  of 
the  cookie-cutter  sensor's  radius  at  each  time  step. 

The  author  also  experimented  with  the  Poisson  course 
change  rate  A,  to  produce  searcher  and  target  motion  most 
closely  satisfying  the  Random  Search  assumptions. 
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V=200,  U=0 


Figure  8 .  iy(Y)  for  Various  Values  of  A 

„  2F 

As  indicated  m  Figure  8,  a  good  value  is  A  -  , — =  . 

Vdd 

This  formula  implies  that  on  the  average,  two  course  change 
events  occur  during  the  time  required  for  the  mover  to  go 

from  edge  to  edge.  If  ^  j_s  too  small,  the  mover  has  a  very 
small  chance  to  change  course  before  bouncing  off  the  edge. 

On  the  other  hand,  if  ^  is  too  large  then  the  mover' s 
position  will  be  potentially  limited  to  a  small  part  of  the 
searcher  and  target  area  (see  Figure  9) . 
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Figure  9.  Searcher's  Movement  Pattern  for  Various 

Values  of  A 


3 .  Output 

Using  the  recommended  values  for  the  reflection 

„  2V 

parameter  (0.5  radian)  and  course  change  rate  (  A  =  - - ==■  ), 

the  Random  Search  simulation  produced  results  very  close  to 
the  expected  theoretical  model . 
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Figure  10.  FT(t)  for  Various  Target  Speeds 
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In  Figure  10,  the  pink  line  is  plots  of  following 
formula : 

Ft (?)  =  1  -  (1  - 7i R2  /  AJx  exp(-2R Vt /  As )  .  (16) 

The  blue  line  represents  simulation  results.  The 
results  indicate  that  this  MATLAB  simulation  is  a  valid 
representation  of  Random  Search. 

It  is  also  clear  that  the  movement  of  the  target 
increases  the  opportunity  for  detection.  Thus,  if  the 
searcher  conducts  a  Random  Search,  the  best  strategy  for  the 
target  is  to  be  stationary. 
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IV.  SIMULATION  OF  THE  EXTENDED  RANDOM  SEARCH  MODEL 
WHERE  SEARCHER  AND  TARGET  AREAS  ARE  NOT  COINCIDENT 


A.  DESCRIPTION  OF  THE  GENERALIZED  RANDOM  SEARCH  MODEL 

1.  Characteristics  of  the  Searcher  and  the  Target 

The  only  difference  from  the  previous  model  is  that  the 
search  area  is  now  not  completely  overlapped  with  the  target 
area . 

B.  COMPUTER  ALGORITHM 

1.  Inputs  [Units] 

•  Number  of  simulation  replications,  N  = 500. 

•  Maximum  simulation  time,  ?max=150  [hour]. 

•  The  length  of  search  area  in  X  direction,  /v  =150 

[nm]  . 

•  The  length  of  search  area  in  Y  direction,  /v  =  1 50 

[nm]  . 

•  The  length  of  target  area  in  X  direction,  =150 

•  [nm]  . 

•  The  length  of  target  area  in  Y  direction,  /v?  =150 

[nm]  . 

•  Searcher  speed,  V =200  [nm/hour] . 

•  Target  speed,  U  [nm/hour] . 

•  Searcher' s  detection  range,  R=  2  [nm]  . 

•  Searcher's  course  change  rate,  As  [1/hour]. 

•  Target's  course  change  rate,  At  [1/hour]  . 
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The  unit  time  of  simulation.  At  = - =  0.005  [hour] 

2F 

The  size  of  search  area,  Ms  =  1502  [  nm2  ]  . 

The  size  of  target  area,  ^=  1502  [  nm2  ]  . 

The  size  of  overlap  area,  ^=1002  [  nm2  ]  . 

Functioning  of  the  Program 


When  this  model  was  implemented  in  MATLAB,  the  faster 
mover's  At  was  used.  For  example,  if  the  searcher's  speed 
is  200kts  and  the  target's  speed  is  lOOkts  then,  this 

2nm 


simulation' s  At  is 


2x  200kts 


=  0.005  hour 


Figure  11 


Target  Movement  Behavior  for  At  (V=200kts, 
U=2kts) 


The  image  in  Figure  11  results  when  using  the  following 
parameters:  tmix  =2  00  hours,  search  area  ^s=  1502nm2  ,  target 


moving  area  J?=1002nm2  ,  overlap  area  Ast  =  502nm2  ,  V  =200kts, 

U  =  2kts  ,  R  =  2nm  ,  At  =  — —  =  0.005  hours  ,  X  =  ^  x  200kts  _  ^  ^  ^ 

2  x  200kts  Vl50nmxl50nm 

,  2  x  2kts 

,  and  A  =  , —  = =0 . 04 /hour .  For  both  the  fast  searcher 

VlOOnmxlOOnm 

and  the  slow  target,  there  is  an  average  of  two  course 
changes  during  the  time  the  mover  crosses  its  area. 


3 .  Output 


Models  for  the  fast  and  slow  target  motion  were 
introduced  in  Chapter  II.  Figure  12  shows  simulation  results 
for  several  target  speeds  plotted  against  the  fast  and  slow 
target  models: 

1.  Green  line  :  fast  target  assumption 

FTfast  (0  =  1 -  exp(-2R  V(Ast  /  A,  )t  /  As )  (17) 

2.  Pink  line  :  slow  target  assumption 

Fmow{t)  =  ( 4  /  4  )(1  -  exp(-2i?  ft/AJ)  (18) 

3.  Blue  line  :  simulation  results 
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Figure  12.  FT(t)  for  Various  Target  Speed,  Overlap  Area 

Asl  =  1 00  nm  x  1 00  nm 


The  fast  target  and  slow  target  equations  can  be 
combined  as  follows  to  approximate  the  simulation  results 
for  any  target  speeds : 


FTcombined  (0  =  (exp (-%4t))FTslm{t)  +  (1  ~  eXP (~^))Frfast  (t) 

a  1  -a 

where  ^  r  [  1  /  Vtime  ]  . 

jAsAtRV 


(19) 
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a  decreases  from  1  to  0  as  the  target  speed  U  increases 
and  as  time  goes  to  infinity,  resulting  in  the  fast  target 
equation  (17)  .  Figure  13  shows  plots  of  simulation  results 
and  FTcombined(t)  .  The  red  line  in  the  Figure  13  represents  the 
results  of  equation  (19). 


Figure  13 


The  Comparison  between  Simulation  Results  and 


F 


Tcombined 


(0 


Given  these  results,  it  appears  that  equation  (19)  can 
be  used  as  a  conservative  estimate  (that  is,  a  lower  bound) 
of  FT(t )  when  searcher  and  target  patrol  areas  are  not 
identical . 
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V.  SIMULATION  OF  A  SYSTEMATIC  RANDOM  SEARCH 
AGAINST  A  RANDOMLY  MOVING  TARGET 


A.  DESCRIPTION  OF  'SYSTEMATIC  SEARCH '  MODEL 

Actual  searchers  cannot  perform  Exhaustive  Search,  but 
it  is  often  approximated  by  parallel  sweeps,  which  is 
similar  to  mowing  the  lawn.  This  is  called  'Systematic 
Search' .  Other  forms  of  Systematic  Search  include  spiral-in 
and  spiral-out  tracks. 


(1)  (2) 

Figure  14.  Simple  Examples  of  'Systematic  Search' 

Patterns 

Figure  14  illustrates  simplified  examples  of  a 
'Systematic  Search'  track.  We  examine  the  two  types  of  going 
and  returning  'Systematic  Search'  patterns.  One  is  a 
'Reverse  Course  Systematic  Search'  in  which  the  searcher 
goes  back  via  the  exact  same  track  to  the  starting  point. 
The  other  is  a  'Cross  Course  Systematic  Search'  in  which  the 
searcher  follows  the  track  shown  in  Figure  14(1)  and  then 

goes  back  via  the  Figure  14(2)  track. 
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B.  COMPUTER  ALGORITHM 

1.  Input  [Units] 

•  Number  of  simulation  replications,  N  =1000. 

•  Maximum  simulation  time,  £max  =54.72  [hour]. 

•  The  length  of  search  area  of  X  direction,  /v  =  148 
[nm]  . 

•  The  length  of  search  area  of  Y  direction,  /v  =  148 
[nm]  . 

•  Searcher  speed,  V =200  [nm/hour] . 

•  Target  speed,  U  [nm/hour] . 

•  Searcher's  detection  range,  R =2  [nm] . 

2  U 

•  Target  course  change  rate,  At  =  — =  [1/hour]. 

V4x/>’ 

•  The  unit  time  of  simulation.  At  = - =  0.005  [hour]  . 

2F 

•  The  searcher' s  moving  distance  during  each  time 
s t 1 ep ,  dunitJime  =  V  x  At  =  1  [  nm  ]  . 

•  The  size  of  search  area,  =  1482  [  nm2  ]  . 

2.  Comparison  of  the  Searcher's  Movement  Patterns 

Simulation  results,  illustrated  in  Figure  15,  indicate 
very  little  difference  between  the  two  search  paths.  The 
author  postulates  that  any  reasonable  Systematic  Search  plan 
(e.g.  parallel  sweeps,  spiral-in,  or  spiral-out)  will  be 
effective  as  long  as  the  searcher  attempts  to  cover  all 
points  in  the  search  area  with  the  sensor  one  time  before 
any  point  is  covered  twice. 
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V=200,  U=2  V=200,  U=2 


Simulation  Result 
Random  Search 
Exhaustive  Search 


Simulation  Result 
Random  Search 
Exhaustive  Search 


V=200,  U=4  V=200,  U=4 


Simulation  Result 
Random  Search 
Exhaustive  Search 


Simulation  Result 
Random  Search 
Exhaustive  Search 


Figure  15.  Comparison  Between  Reverse  and  Cross  Course 


In  Figure  15,  each  line  is  defined  as  follows: 
Exhaustive  Search  (Black  line)  :  FT(f)  =  (2VRt)/As  . 

Random  Search  (Pink  line)  :  FT(t)  =  X-il-nR1  / As)xexp(-2RVt / As)  . 
Simulation  Results (Blue  line) . 

As  expected,  there  is  very  little  difference  between 
the  reverse  and  the  cross  course's  simulation  results. 
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3 .  Output 
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Figure  16. 


FT{f) 


for  Various  Target  Speed 


(2  Times  Sweep) 
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In  general,  it  is  thought  that  a  simulation  of 

Systematic  Search  always  performs  better  than  Random  Search. 
Figure  16  suggests  that  this  is  true. 

The  measure  of  effectiveness  used  when  evaluating 

search  plans  was  "Incomplete  Mean  Time  to  Detection 
evaluated  at  time  t  (  IMTDO )  )  ,  "  defined  as , 

IMTD(T)  =  f'(l  -FT(x))dx.  IMTD(0  is  used  as  a  surrogate  for 

J  0 

Mean  time  to  detection  (  E(T )  )  .  In  fact,  limIMTD(d )  =  E(T).  We 

t  ->oo 

use  IMTD(7  )  because  we  often  cannot  run  the  search 

simulation  long  enough  to  accurately  estimate  E(T )  (see 
Figure  17 ) . 


09 

08 

07 

06 
0.5 
04 

03 

03 

01 

0 

I 

Time 

Figure  17.  Example  of  'Incomplete  Mean  Time  to 

Detection'  where  t  =50  hours. 


V=200.  U=20 


Area  1 


Area  2 


Area  1  +  Area  2  =  f  (1-  Fr{.\))<h  =  E(T)  -  Mean  tune  to  detection 
JitiO  1 

iso 

^  (1-  Fr(y))dx  =  Incomplete  mean  tune  to  detection 
Incomplete  mean  tune  to  detection  <  Mean  tune  to  detection 
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We  assume  that  the  target  would  like  to  select  a  speed 
maximizing  IMTD(f)  .  In  Figure  18,  IMTD(0  is  dramatically 
increased  when  the  target  speed  changed  from  Okts  to  8kts. 
And  there  are  no  significant  changes  between  8kts  and  44kts. 
For  target  speeds  greater  than  about  44kts,  IMTD(d)  steadily 
decreases . 


Systematic  Search(V=200kts.  A=148nm'148nm) 


Figure  18.  Incomplete  Mean  Time  to  Detection  for  Various 

T 

Naval  Postgraduate  School  Professor,  Alan  R.  Washburn, 
suggested  this  when  he  wrote,  "a  target  that  wishes  to  avoid 
detection  might  actually  choose  to  move  around  atU  =  .2V,  on 
the  grounds  that  this  is  enough  motion  to  prevent  an 
Exhaustive  Search,  but  nonetheless  increases  the  equivalent 
searcher  speed  by  only  1%."2  After  that  speed,  as  Professor 
Washburn  expected,  the  IMTD(7)  decreases  because  the  fast 


2  Alan  R.  Washburn,  Search  and  detection,  4th  ed.  6-3. 
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random  movement  of  the  target  runs  it  into  the  searcher  more 
often  than  away  from  the  searcher.  Figure  18  shows  that  even 
though  the  speed  of  target  increases  from  Okts  to  40kts, 
which  is  a  speed  enough  to  significantly  degrade  an 

Exhaustive  Search,  the  increase  of  V  is  only  1%  of  searcher 
speed . 

All  things  considered,  the  author  suggests  that  the 
optimal  evasion  speed  for  target  evading  a  Systematic  Search 
is  approximately  0.05  to  0.2  of  the  searcher  speed. 


V=200kts 


Figure  19.  The  Relationship  Between  Target  Speed ( U  )  and 

V 
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VI .  CONCLUSIONS 


A.  CONCLUSIONS  AND  RECOMMENDATIONS 

The  main  contributions  of  this  thesis  can  be  narrowed 
down  to  two  results.  The  first  result  is  a  generalization  of 
the  Random  Search  formula  which  allows  both  searcher  and 
target  to  move  and  does  not  require  the  patrol  areas  to  be 
identical . 

The  other  result  is  that  the  best  speed  for  a  randomly 
moving  target  evading  a  Systematic  Search  ranges  from  0.05 
to  0.2  of  the  searcher  speed. 

In  this  thesis,  the  searcher  always  uses  a  cookie- 
cutter  sensor.  We  could  develop  a  potentially  more  realistic 
model  by  relaxing  the  cookie-cutter  sensor  assumption  in 
future  studies.  For  instance,  if  we  assume  the  probability 
of  target  detection  with  sensor  range  R  is  0.8,  then 

f  1,  if  Uniform  _  Random (0, 1)  <  0.8  &  range  <  R 
.P(detection)=P/  =  <  . 

[0,  otherwise  . 

In  addition,  the  comparison  of  results  when  using  various 
kinds  of  sensor  types  such  as  cookie-cutter  sensor,  inverse- 
cube  law  sensor,  and  triangle  sensor  would  be  worthwhile  to 
study . 
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APPENDIX 


A.  MATLAB  CODE 


1 .  Random  Search  Model 


Nreps=500 ; 
tmax=150; 
lx=150 ; 
ly=150; 

V=200; 

U=50  ; 

R=2  ; 

dt=R/ (2*max (V, U) ) ; 
lams= (2*V) / sqrt (lx*ly)  ; 
lamt= (2*U) / sqrt (lx*ly)  ; 
Xs= zeros ( 1 , tmax/dt+1 ) ; 
Ys=Xs ; 

Cs=Xs ; 

Xt=Xs ; 

Yt=Xs ; 

Ct=Xs ; 

CumDet=Xs ; 

T=0 : dt : tmax; 

A=lx*ly; 

for  n=l : Nreps 

xs=rand* (lx-2*R) +R; 

ys=rand* (ly-2*R) +R; 

xt=rand*lx; 

yt=rand*ly; 


%number  of  simulation  replications 
%max .  simulation  time  (hr) 

%search  area  length  in  x  direction  (nm) 
%search  area  length  in  y  direction  (nm) 
%searcher  speed  (nm/hr) 

%target  speed  (nm/hr) 

%searcher  detection  range  (nm) 

%delta  t  (hours) 

%searcher  course  change  rate  (1/hr) 

%target  course  change  rate  (1/hr) 

%initialize  x-position  to  zero ( searcher ) 
%initialize  y-position  to  zero ( searcher ) 
%initialize  searcher  course  to  zero 
%initialize  x-position  to  zero (target) 
%initialize  y-position  to  zero (target) 
%initialize  target  course  to  zero 
%initialize  cumulative  detection  stats 
%simulation  time  vector 
%search  area 
%main  simulation  loop 

%initial  searcher  and  target  x  and  y  positions 


cs=rand*2*pi; 

ct=rand*2*pi; 

t=0; 

tindex=l ; 

Xs (tindex) =xs; 

Ys (tindex) =ys; 

Cs (tindex) =cs; 

Xt (tindex) =xt; 

Yt (tindex) =yt; 

Ct (tindex) =ct; 
for  t=l:tmax/dt 

tindex  =  tindex+1; 
if  rand<lams*dt; 
cs=rand*2  *pi ; 


initial  searcher  course 
initial  target  course 
set  simulation  time  to  0 
initialize  time  index  to  1 
save  initial  searcher  x  position 
searcher  y  position 
searcher  course 
target  x  position 
target  y  position 
target  course 


%save 

%save 

%save 

%save 

%save 

%inner 


initial 

initial 

initial 

initial 

initial 

loop 


%update  simulation  time  index 
%Poisson  course  change  rate. 
%New  courses  uniform  (0,2*pi) . 


end 

if  Xs (tindex-1 ) <R;  cs= ( rand- . 5 ) ; end 
if  Xs (tindex-1) > (lx-R) ;  cs=pi+ (rand- . 5) ; end 
if  Ys (tindex-1 ) <R;  cs=pi/2+ (rand- . 5) ; end 
if  Ys (tindex-1) > (ly-R) ;  cs=-pi/2+ (rand- . 5) ; end 

Xs (tindex)  =  Xs (tindex-1) +V*dt*cos (cs) ;  %Update  x  and  y  positions 
Ys (tindex)  =  Ys (tindex-1) +V*dt*sin (cs) ; 
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Cs (tindex) =cs; 

if  rand<lamt*dt;  ct=rand*2 *pi ; end 
if  Xt (tindex-1) <0;  ct= ( rand- . 5 ) ; end 
if  Xt (tindex-1) >lx;  ct=pi+ (rand- . 5) ; end 
if  Yt (tindex-1) <0;  ct=pi/2+ (rand- . 5) ; end 
if  Yt (tindex-1) >ly;  ct=-pi/2+ ( rand- . 5 ) ; end 

Xt (tindex)  =  Xt (tindex-1 ) +U*dt*cos  (ct) ;  %Update  x  and  y  positions 
Yt (tindex)  =  Yt (tindex-1 ) +U*dt*sin (ct) ; 

Ct (tindex) =ct; 

end  dinner  loop  (time  increasing  from  0  to  tmax) 

CumDet  =  CumDet  +  cummax ( (Xs-Xt) .  A2  +  (Ys-Yt) .A2  <=  RA2); 
end  %outer  loop  (simulation  replications) 

Probability=CumDet . /Nreps ; 

RSfun=l- (l-pi*R*R/A) *exp ( -2 *R*vtilde (V, U) *T/A) ; 

plot (T, CumDet/Nreps ,  ' b- '  ,  T,  RSfun,  ' m- ' ) ,  axis ([ 0 , tmax, 0 , 1 ] ) 
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2. 


Patrol  Areas  are  not  Identical'  Model 


Nreps=500 ; 

tmax=150; 

lx=150 ; 

ly=150; 

lx_t=150; 

ly_t=150; 

V=200; 

U=  1 5 ; 

R=2  ; 

dt=R/ (2*max (V, U) ) ; 
lams= (2*V) / sqrt (lx*ly)  ; 
lamt= (2*U) / sqrt (lx_t*ly_t) ; 
A=lx*ly; 

Ast=100*100; 

At=lx_t*ly_t; 

Xs= zeros ( 1 , tmax/dt+1 ) ; 


%search  area  length  in  x  direction  (nm) 

%search  area  length  in  y  direction  (nm) 

%target  area  length  in  x  direction  (nm) 

%target  area  length  in  y  direction  (nm) 

%searcher  speed  (nm/hr) 

%target  speed  (nm/hr) 

%searcher  detection  range  (nm) 

%delta  t  (hours) 

%searcher  course  change  rate  (1/hr) 
%searcher  course  change  rate  (1/hr) 
%search  area (As) 

%overlap  area 
%target  area (At) 


Ys=Xs ; 

Cs=Xs ; 

Xt=Xs ; 

Yt=Xs ; 

Ct=Xs ; 

CumDet=Xs ; 

T=0 : dt : tmax; 
for  n=l : Nreps 
xs=rand* (lx-2*R) +R; 
ys=rand* (ly-2*R) +R+50; 
xt=rand*lx  t+50; 
yt=rand*ly_t; 
cs=rand*2  *pi ; 
ct=rand*2*pi; 
t=0; 

tindex=l ; 


Xs (tindex) =xs; 

%save 

initial 

searcher 

x  position 

Ys (tindex) =ys; 

%save 

initial 

searcher 

y  position 

Cs (tindex) =cs; 

%save 

initial 

searcher 

course 

Xt (tindex) =xt; 

%save 

initial 

target  x 

position 

Yt (tindex) =yt; 

%save 

initial 

target  y 

position 

Ct (tindex) =ct; 

%save 

initial 

target  course 

for  t=l:tmax/dt 

tindex  =  tindex+1; 
if  rand<lams*dt; 
cs=rand*2*pi; 
end 

if  Xs (tindex-1 ) <R;  cs= (rand- . 5) ; end 
if  Xs (tindex-1) > (lx-R) ;  cs=pi+ (rand- . 5) ; end 
if  Ys (tindex-1) < (R+50) ;  cs=pi/2+ (rand- . 5) ; end 
if  Ys (tindex-1) > (ly-R+50) ;  cs=-pi/2+ (rand- . 5 ) ; end 
Xs (tindex)  =  Xs (tindex-1) +V*dt*cos (cs) ; 

Ys (tindex)  =  Ys (tindex-1) +V*dt*sin (cs) ; 

Cs  (tindex) =cs; 

if  rand<lamt*dt;  ct=rand*2 *pi ; end 

if  Xt (tindex-1) <50;  ct= (rand- . 5) ; end 

if  Xt (tindex-1) > (lx  t+50);  ct=pi+ ( rand- . 5 ) ; end 
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if  Yt (tindex-1) <0;  ct=pi/2+ (rand- . 5) ; end 
if  Yt (tindex-1) >ly  t;  ct=-pi/2+ (rand- . 5 ) ; end 
Xt(tindex)  =  Xt (tindex-1 ) +U*dt*cos (ct) ; 

Yt(tindex)  =  Yt (tindex-1 ) +U*dt*sin (ct) ; 

Ct  (tindex) =ct; 
end 

CumDet  =  CumDet  +  cummax ( (Xs-Xt) .  A2  +  (Ys-Yt) .A2  <=  RA2); 
if  n/50  ==  floor(n/50),  n,  end; 

end 

RSfun_slow= (Ast/At) * (1-exp (- (2*R*vtilde (V, U) *T) /A) ) ; 
RSfun_fast=l-exp (- (2*R* (Ast/At) *vtilde (V,U) *T) /A) ; 
subplot (2 , 1 , 1 ) 

plot (T, CumDet/Nreps , ' b- ' ,  T,RSfun  slow, ' m- ' , T, RSfun  fast, 'g'), 
axis ( [ 0 , tmax, 0,1]) 
subplot  (2 , 1,2) 

xline= [150  0  0  150  150  50  50  200  200  50  50]; 
yline= [50  50  200  200  50  50  150  150  0  0  50]; 
plot(Xs,Ys,  'b-',Xt,Yt,  'r 

'  ,  Xs ( 1 )  ,  Ys ( 1 )  ,  ' yo ' , Xt ( 1 ) , Yt ( 1 ) ,  ' yo ' , xline, yline,  ' k ' ) , axis  off 
hold  on 

axis([0  lx+50  0  ly+50] ) ,  axis  square 
scale=lx+50 ; 

text (Xs(l)+scale/40,Ys(l), 'Start', ' Font Size ' , 10, ' FontName ' , ' Times ' ) 
text (Xt(l)+scale/40,Yt(l),  'Start',  ' Font Size  '  ,  10,  ' FontName ' ,  ' Times ' ) 
hold  off 
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3. 


Systematic  Search  Model 


Nreps=l 000 ; 

lx=148; 

ly=148; 

V=200; 

U=150; 

R=2  ; 

dt=R/ (2*V) ; 

unit  time  movement=V*dt; 
tmax=dt  *10944; 
lams= (2*V) / sqrt (lx*ly) ; 
lamt= (2*U) / sqrt (lx*ly) ; 
Xs=zeros ( 1 , tmax/dt+1 )  ; 
Ys=Xs ; 

Cs=Xs ; 

Xt=Xs ; 

Yt=Xs ; 

Ct=Xs ; 

CumDet=Xs ; 

CumDet  s=Xs; 

T=0 : dt : tmax; 

A=lx*ly; 

O, _ 

o 

xs_s=R; 
for  h=l : 36 


%number  of  simulation  replications 
%search  area  length  in  x  direction  (nm) 
%search  area  length  in  y  direction  (nm) 
%searcher  speed  (nm/hr) 

%target  speed  (nm/hr) 

^searcher  detection  range  (nm) 

%delta  t  (hours) 

%max .  simulation  time  (hr) =54. 72  hours 
%searcher  course  change  rate  (1/hr) 
%target  course  change  rate  (1/hr) 
%initialize  x-position  to  zero ( searcher ) 
%initialize  y-position  to  zero ( searcher ) 
%initialize  searcher  course  to  zero 
%initialize  x-position  to  zero (target) 
%initialize  y-position  to  zero (target) 
%initialize  target  course  to  zero 
%initialize  cumulative  detection  stats 
%initialize  cumulative  detection  stats 
%simulation  time  vector 
%search  area 

define  Searcher' s  coordinate 
%initial  x  coordinate ( searcher ) 


xs  column=linspace (4*h-2, 4*h-2, 144)  ;  xs  s=[xs  s  xs  column] ; 
xs  row= (4*h-l) :unit  time  movement :( (2*R) /unit  time  movement+4 *h-2 ) 
xs  s=[xs  s  xs  row]; 

end 


xs  f=linspace (146,  146,  144)  ; 
xs_s=[xs_s  xs_f ] ; 
for  b=36 : -1 : 1 

xs  column  b=linspace (4*b+2, 4*b+2, 144) ;  xs  s=[xs  s  xs  column  b] ; 
xs  row  b=(4*b+l) :-unit  time  movement: (4*b+2- 
(2*R)/unit  time  movement); 
xs  s=[xs  s  xs  row  b] ; 

end 

xs_b=linspace (2,2,144)  ; 
xs_s=[xs_s  xs_b]; 

ys_s=2 : 1 : 146;  %define  target's  coordinate 

for  z=l : 1 8 

ys_row=linspace (146, 146, 4) ;  ys_s=[ys_s  ys_row] ; 
ys_column=145 : -1 : 2 ;  ys_s=[ys_s  ys_column]; 
ys_row=linspace (2 , 2 , 4 ) ;  ys_s=[ys_s  ys_row] ; 
ys_column=3 : 1 : 146;  ys_s=[ys_s  ys_column] ; 
end 

for  z=l:18 

ys_column=145 : -1 : 2 ;  ys_s=[ys_s  ys_column]; 
ys_row=linspace (2 , 2 , 4 ) ;  ys_s=[ys_s  ys_row] ; 
ys_column=3 : 1 : 146;  ys_s=[ys_s  ys_column]; 
ys_row=linspace (146,  146,  4) ;  ys_s=[ys_s  ys_row] ; 
end 

ys_b=145 : -1 : 2 ;  ys_s=[ys_s  ys_b] ; 

o, _ 

o 
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for  n=l : Nreps 
xs=rand* (lx-2*R) +R; 
ys=rand* (ly-2*R) +R; 
xt=rand*lx; 
yt=rand*ly; 
cs=rand*2  *pi ; 
ct=rand*2  *pi ; 
t=0; 

tindex=l ; 

Xs (tindex) =xs; 

Ys (tindex) =ys; 

Cs (tindex) =cs; 

Xt (tindex) =xt; 

Yt (tindex) =yt; 

Ct (tindex) =ct; 
for  t=l:tmax/dt 


%main  simulation  loop 

%initial  searcher  and  target  x  and  y  positions 


tindex  =  tindex+1; 
if  rand<lams*dt; 
cs=rand*2  *pi ; 


%initial  searcher  course 
%initial  target  course 
%set  simulation  time  to  0 
%initialize  time  index  to  1 
%save  initial  searcher  x  position 
%save  initial  searcher  y  position 
%save  initial  searcher  course 
%save  initial  target  x  position 
%save  initial  target  y  position 
%save  initial  target  course 
%inner  loop 

%update  simulation  time  index 
%Poisson  course  change  rate. 

%New  courses  uniform  (0,2*pi) . 


end 

if  Xs (tindex-1 ) <R;  cs= ( rand- . 5 ) ; end 
if  Xs (tindex-1) > (lx-R) ;  cs=pi+ (rand- . 5) ;end 
if  Ys (tindex-1) <R;  cs=pi/2+ (rand- . 5 ) ; end 
if  Ys (tindex-1) > (ly-R) ;  cs=-pi/2+ (rand- . 5) ; end 
Xs (tindex)  =  Xs (tindex-1) +V*dt*cos (cs) ; 

Ys (tindex)  =  Ys (tindex-1) +V*dt*sin (cs) ; 

Cs  (tindex) =cs; 

if  rand<lamt*dt;  ct=rand*2 *pi ; end 
if  Xt (tindex-1) <0;  ct= (rand- . 5) ; end 
if  Xt (tindex-1) >lx;  ct=pi+ ( rand- . 5 ) ; end 
if  Yt (tindex-1) <0;  ct=pi/2+ (rand- . 5) ; end 
if  Yt (tindex-1) >ly;  ct=-pi/2+ (rand- . 5 ) ; end 
Xt (tindex)  =  Xt (tindex-1 ) +U*dt*cos (ct) ; 

Yt (tindex)  =  Yt (tindex-1 ) +U*dt*sin (ct) ; 

Ct (tindex) =ct; 

end  %inner  loop  (time  increasing  from  0  to  tmax) 

CumDet  =  CumDet  +  cummax ( (Xs-Xt) .A2  +  (Ys-Yt) .A2  <=  RA2); 

CumDet  s  =  CumDet  s  +  cummax ((xs  s-Xt) .A2  +  (ys  s-Yt) . A2  <=  RA2); 
end  %outer  loop  (simulation  replications) 

Exhautive  Search=2 *R*vtilde (V, U) *T/A; 
subplot  (2 , 1 , 1 ) 

plot (T, CumDet/Nreps, 'b-',  T, CumDet_s/Nreps, ' m- 
' , T, Exhautive_Search,  ' k'  )  ,  axis ( [0, tmax, 0,1]) 
subplot (2 , 1,2) 

plot (xs_s , ys_s , Xt , Yt , ' r- ' , xs_s ( 1 ) , ys_s ( 1 ) , ' ko ' , Xt ( 1 ) , Yt ( 1 ) , ' go ' ) , axis ( [ 0 
lx  0  ly]),axis  square 
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